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Abstract: 

The radiative corrections of the strong and electroweak interactions are calculated 
at next-to-leading order for Higgs-boson production in the weak-boson-fusion channel at 
hadron colliders. Specifically, the calculation includes all weak-boson fusion and quark- 
antiquark annihilation diagrams to Higgs-boson production in association with two hard 
jets, including all corresponding interferences. The results on the QCD corrections confirm 
that previously made approximations of neglecting s-channel diagrams and interferences 
are well suited for predictions of Higgs production with dedicated vector-boson fusion 
cuts at the LHC. The electroweak corrections, which also include real corrections from 
incoming photons and leading heavy-Higgs-boson effects at two-loop order, are of the same 
size as the QCD corrections, viz. typically at the level of 5—10% for a Higgs-boson mass 
up to ~ 700 GeV. In general, both types of corrections do not simply rescale differential 
distributions, but induce distortions at the level of 10%. The discussed corrections have 
been implemented in a flexible Monte Carlo event generator. 
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1 Introduction 



The production of a Standard Model Higgs boson in association with two hard jets in 
the forward and backward regions of the detector — frequently quoted as the "vector-boson 
fusion" (VBF) channel — is a cornerstone in the Higgs search both in the ATLAS [1] and 
CMS [2] experiments at the LHC. This is not only true for the Higgs-mass range between 
100 and 200 GeV, which is favoured by the global Standard Model fit to electroweak (EW) 
precision data [3], but also for a Higgs mass of the order of several 100 GeV up to the 
theoretical upper limit set by unitarity and triviality constraints. Higgs production in the 
VBF channel also plays an important role in the determination of Higgs couplings at the 
LHC (see e.g. Ref. [4]). Even bounds on non-standard couplings between Higgs and EW 
gauge bosons can be imposed from precision studies in this channel [5] . 

The production of Higgs+2jets receives two kinds of contributions at hadron colliders. 
The first type, where the Higgs boson couples to a weak boson that links two quark lines, 
is dominated by squared t- and ti-channel-like diagrams and represents the genuine VBF 
channel. The hard jet pairs have a strong tendency to be forward-backward directed 
in contrast to other jet production mechanisms, offering a good background suppres- 
sion (transverse-momentum and rapidity cuts on jets, jet rapidity gap, central-jet veto, 
etc.). Applying appropriate event selection criteria (see e.g. Refs. [6-10] and references in 
Refs. [11,12]) it is possible to sufficiently suppress background and to enhance the VBF 
channel over the second Higgs-|-2jets production mechanism that mainly proceeds via 
strong interactions. In this second channel the Higgs boson is radiated off a heavy-quark 
loop that couples to any parton of the incoming hadrons via gluons [13,14]. According to 
a recent estimate [15] hadronic production contributes about 4—5% to the Higgs+2jets 
events for a Higgs mass of 120 GeV after applying VBF cuts. A next-to-leading order 
(NLO) analysis of this contribution [14] shows that its residual scale dependence is still 
of the order of 35%. 

Higgs production in the VBF channel is a pure EW process in leading order (LO) 
involving only quark and antiquark parton distributions. As s-channel diagrams and 
interferences tend to be suppressed, especially when imposing VBF cuts, the cross section 
can be approximated by the contribution of squared t- and w-channel diagrams only. The 
corresponding QCD corrections reduce to vertex corrections to the weak-boson-quark 
coupling. Explicit NLO QCD calculations in this approximation [11,16-19] confirm the 
expectation that these QCD corrections are small, because they are shifted to the parton 
distribution functions (PDFs) via QCD factorization to a large extent. The resulting 
QCD corrections are of the order of 5—10% and reduce the remaining factorization and 
renormalization scale dependence of the NLO cross section to a few per cent. 

In a recent letter [20] we completed the existing NLO calculations for the VBF channel 
in two respects. Firstly, we added the full NLO EW corrections. Secondly, we calculated 
the NLO QCD corrections including, for the first time, the complete set of QCD diagrams, 
namely the t-, u-, and s-channel contributions, as well as all interferences. Focussing on 
the integrated cross section (with and without dedicated VBF selection cuts) , we discussed 
the impact of EW and QCD corrections in the favoured Higgs-mass range between 100 
and 200 GeV. We found that the previously unknown NLO EW corrections are of the 
order of —5% and, thus, as important as the QCD corrections. In the EW corrections we 
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Figure 1: Topologies for t-, m-, and s-channel contributions to qq ggH in LO, where q 
denotes any quark or antiquark and V stands for W and Z bosons. 

also take into account real corrections induced by photons in the initial state and QED 
corrections implicitly contained in the DGLAP evolution of PDFs. We found that these 
photon-induced processes lead to corrections at the per-cent level. 

In this paper we describe more details of our calculation, which is performed in a widely 
analogous way to the EW and QCD corrections to the Higgs decay H WW/ZZ 
4fermions [21,22]. We classify the NLO QCD corrections into four different categories; 
the previously known corrections [11,16-19] are contained in one of these categories. More- 
over, we extend our numerical discussion of the EW and QCD corrections in two respects. 
We now consider cross sections for Higgs masses above 200 GeV, including the leading EW 
two- loop corrections oc G^M^ for a heavy Higgs boson using the results of Ref. [23], and 
we discuss differential distributions in transverse momenta, in rapidities, in the azimuthal 
angle difference of the tagging jets, and in the jet-jet invariant mass. We pay particu- 
lar attention to the issue of distortions in distributions induced by radiative corrections, 
because such distortions usually are the signature of non-standard couplings. 

The paper is organized as follows: In Section 2 we describe the calculational setup 
and classify the QCD corrections into different categories. The discussion of numerical 
results is presented in Section 3. Section 4 contains our conclusions. 

2 Details of the calculation 

2.1 General setup 

At LO, the hadronic production of Higgs+2jets via weak bosons receives contributions 
from the partonic processes qq — > ggH, qq ggH, qq qqH, and qq qqR. For each 
relevant configuration of external quark flavours one or two of the topologies shown in 
Figure 1 contribute. All LO and one-loop NLO diagrams are related by crossing symmetry 
to the corresponding decay amplitude H — > qqqq. The QCD and EW NLO corrections to 
these decays were discussed in detail in Refs. [21,22], in particular a representative set of 
Feynman diagrams can be found there. 

To be more specific, we first show how the lowest-order and loop amplitudes for sub- 
processes of the type qq qqB. can be obtained from the corresponding results for 
H qaQbqcqd- The basic lowest-order decay amplitudes involving W- or Z-boson ex- 
change, called Mo^''''''''''^''^{ka, h, K, ka) with V = W, Z, have been defined in Eq. (2.8) 
of Ref. [21]; the two potentially relevant tree diagrams are shown in Figure 2; the cor- 
responding squares and interference are illustrated in Figure 3. The external momenta 
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Figure 2: Generic lowest-order diagrams for the Higgs decay H —>■ qaQbClcQd, where V, V = 
W, Z denote the exchanged weak bosons. The lowest-order diagrams for qq qqB. are 
obtained by crossing any pair of (anti-)quarks into the initial state and the Higgs boson 
into the final state. 




Figure 3: Types of squared and interference diagrams contributing in lowest order. 

{pi,p'i}, helicities {a^, a^}, and colour indices {q, c-} are assigned to the scattering particles 
according to 

q{pu c^i, ci) + q{p2, C2) q{p[, a[, c[) + q{p2, a'^, 4) + B.{p'^). (2.1) 

In order to compactify notation, we omit the labels pi, ai, q, etc. in the amplitudes ^^g^w 
of the scattering process, i.e. we imphcitly have = •^Tic^l[i^^"^'^^'^^{Pi^P2,Pi,P2)^ 

and abbreviate the helicity and momentum assignment in JH^^ as A1^^(122'l') = 
^yy,-cri,-(T2,cr2,o-i(^_p^^ _p2,p2,Pi), etc. Note that momenta and helicities crossed into 
the initial state receive a sign change. In this notation the lowest-order amplitudes ^0 
for the six basic flavour channels in qq qqB. read 

^u,d,^u,d, ^ v^^^.v^;^c'r'^^'=^'^'^A^^W(122lO (11'2'2), 

^u.u,^d,d, ^ s^j6MCT''''''''Mfil22'l') - V,kV;Ct'''''''''M^'^{n'2'2), 

^id,^u,u, ^ 5^^5^^ct'-'^--"'Mf{2'l'l2) - V,*V5,Ci^^^'=^'='^A<^W(2'2ii'), 

jQu,^u,u, ^ 5^^.5^^c'r'=^'^'='^A^o'^(122'lO -5.fc(5,7C'r''^'^'^^o^^(ll'2'2), 

^Jd,^4<i, ^ 5ij5kiCl'''"''"''Mf {122'!') - dik^jiCl'"''"''"" Mf{ll'2'2), (2.2) 
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where i,j,k,l are generation indices, Vij are quark-mixing matrix elements, and C'"'"^'^ is 
one of the two colour operators 




which are relevant to span a general H — > QaqbQcQd amplitude in colour space. The second 
operator Cg^'^'^, which involves the Gell-Mann matrices A'^, becomes relevant in the QCD 
corrections discussed below. The relative sign between the two amplitude contributions 
on the r.h.s. of (2.2) originates from their different fermion- number flow. In Section 2 of 
Ref. [21] the calculation of is described in terms of Weyl-van-der-Waerden spinor 

products {ab) in the conventions of Ref. [24] , where a and b are spinors corresponding to 
external momenta. We note that complex conjugate products (ab)* (but not (ab)) receive 
an additional sign factor for each crossed momentum —pi involved in the product. 

The lowest-order and loop amplitudes for subprocesses of the type qq qqB. and 
qq — > qqB. can be obtained as follows. We assign the external momenta, helicities, and 
colour indices as 

q{pu (Ji, ci) + q{p2, (T2, C2) q{Pi, (t[, q) + qip'^, a'^, 4) + ii{p'^), 

q{pi, ai, Ci) + q{p2, 02, C2) q{p[, a[, c[) + q{p'^, a'^, 4) + ii{p'^). (2.4) 

Then the corresponding amplitudes can be obtained via crossing symmetry from those 
for the process (2.1) as 

,qq^qq,a-ia2a-ia' , / / .qq^qq,-a'.,a2,-(Ti,a', / / 

^qq^qq,'7ia2a[a'2 , ^' ^' \ _ Am^m,<^l-'^2,'^'i-'^2( _ / / . (r, 

When calculating the corresponding cross sections, symmetry factors 1/2 must be taken 
into account for identical fermions or antifermions in the final state. 

In our calculation wc neglect external quark masses whenever possible, i.e. everywhere 
but in the mass-singular logarithms. In (2.2) we made the CKM matrix elements Vij ex- 
phcit. Note that only absolute values of the CKM matrix elements survive after squaring 
the amplitudes; for the squared W-mediated diagrams this is obvious, for the interference 
between W- and Z-mediated diagrams results after contraction of the CKM matrix 

elements with Kronecker deltas. Numerically, only the mixing among the first two gen- 
erations could be relevant, but its impact on Higgs production via VBF was found to be 
negligible. Since the contributions of external b quarks are suppressed, either by bottom 
densities or by s-channel suppression, we optionally include b quarks in the initial and 
final states in our LO predictions, but not in the calculation of corrections. 

2.2 Evaluation of NLO corrections 

Evaluating 2 — > 3 particle processes at the NLO level is non-trivial, both in the 
analytical and numerical parts of the calculation. In order to ensure the correctness of our 
results we have evaluated each ingredient twice, resulting in two completely independent 
computer codes yielding results in mutual agreement. The actual calculation of virtual 
and real NLO corrections for the partonic processes is performed along the same lines as 
described in Ref. [21,22] for the decays H — > 4/. Therefore, we only repeat the salient 
features of the evaluation. 
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(i) Virtual corrections 



The virtual corrections modify the partonic processes that are already present at LO; 
there are about 200 one-loop diagrams per tree diagram in each flavour channel. At NLO 
these corrections are induced by self-energy, vertex, box (4-point), and pentagon (5-point) 
diagrams. The calculation of the EW one-loop diagrams has been performed both in the 
conventional 't Hooft-Feynman gauge and in the background-field formalism using the 
conventions of Refs. [25] and [26], respectively. The QCD one- loop diagrams are evaluated 
in 't Hooft-Feynman gauge. 

In contrast to the t- and w-channel contributions (first two diagrams in Figure 1), the 
s-channel diagrams (last diagram in Figure 1) contain resonant W- or Z-boson propagators 
that require a proper inclusion of the finite gauge-boson widths. For the implementation 
of the finite widths we use the complex-mass scheme, which was introduced in Ref. [27] 
for lowest-order calculations and generalized to the one- loop level in Ref. [28]. In this 
approach the W- and Z-boson masses are consistently considered as complex quantities, 
defined as the locations of the propagator poles in the complex plane. This leads to 
complex couplings and, in particular, a complex weak mixing angle. The scheme fully 
respects all relations that follow from gauge invariance. A brief description of this scheme 
can also be found in Ref. [29]. 

The amphtudes have been generated with FeynArts, using the two independent 
versions 1 and 3, as described in Refs. [30] and [31], respectively. The algebraic evaluation 
has been performed in two completely independent ways. One calculation is based on an 
in-house program written in Mathcmatica, the other has been completed with the help of 
FormCalc [32]. The amplitudes are expressed in terms of standard matrix elements and 
coefficients, which contain the tensor integrals, as described in the appendix of Ref. [33]. 

The tensor integrals are evaluated as in the calculation of the corrections to e+e~ — > 4/ 
[28,34]. They are recursively reduced to master integrals at the numerical level. The 
scalar master integrals are evaluated for complex masses using the methods and results 
of Ref. [35]. UV divergences are regulated dimensionally and IR divergences with an 
infinitesimal photon or gluon mass. Tensor and scalar 5-point functions are directly 
expressed in terms of 4-point integrals [36,37]. Tensor 4-point and 3-point integrals are 
reduced to scalar integrals with the Passarino-Veltman algorithm [38] as long as no small 
Gram determinant appears in the reduction. If small Gram determinants occur, we expand 
the tensor coefficients about the limit of vanishing Gram determinants and possibly other 
kinematical determinants, as described in Ref. [37] in detail. 

Since corrections due to Higgs-boson self-interactions become important for large 
Higgs-boson masses, we have included the dominant two-loop corrections to the l^KH ver- 
tex proportional to G^M^ in the large-Higgs-mass limit which were calculated in Ref. [23]. 
Specifically, we include this effect via a correction factor 



to the squares of the basic LO amplitudes Aig in the t- and w-channel. We do not 
include this correction in the (suppressed) s-channel contributions, because the underlying 
assumption in the derivation of Sq2m4 that Mh is much larger than any other relevant 




(2.6) 
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scale is spoiled by the invariant s that can be of the order of or larger. We do not apply 
^g2m4 to interferences either, because this would require a more complicated structure in 
the correction (involving more than one form factor). The impact of 0{G^M^) corrections 
on interferences and s-channel contributions is certainly negligible, since these effects are 
suppressed themselves. 

(ii) Real corrections 

The matrix elements for the real corrections (photonic/gluonic bremsstrahlung and 
photon-/gluon- induced processes) are obtained via crossing from the bremsstrahlung cor- 
rections to the related Higgs decays, H — > 4/ + 7/g. Explicit amplitudes for H — > 4/ + 7 
are given in Section 4.1 of Ref. [21] in terms of spinor products; for H — 4/ + g such 
results can be found in Section 3.3 of Ref. [22]. The matrix elements relevant for the 
calculation presented here have been checked against results obtained with Madgraph 
[39]. 

The bremsstrahlung corrections involve singularities from soft or collincar pho- 
ton/gluon emission; the photon-/gluon- induced processes contain singularities from 
collincar initial-state splittings. Soft singularities, which arc regularized by an infinitesi- 
mal photon/gluon mass, cancel between virtual and bremsstrahlung corrections. CoUinear 
singularities connected to the initial or final state are regularized by small quark masses, 
which appear only in logarithms. While singularities connected to coUinear configura- 
tions in the final state cancel for "collinear-safe" observablcs automatically after applying 
a jet algorithm, singularities connected to coUinear initial-state splittings are removed via 
factorization by PDF redefinitions, as described in more detail in Section 2.5. 

Technically, the soft and coUinear singularities for real photon emission are isolated 
both in the dipole subtraction method following Ref. [40] and in the phase-space slicing 
method. For photons in the initial state the subtraction and slicing variants described 
in Ref. [41] arc applied. The results presented in the following are obtained with the 
subtraction method, which numerically performs better. 

The phase-space integration is performed with Monte Carlo techniques. One of our two 
codes employs a multi-channel Monte Carlo generator [42] similar to the one implemented 
in Racoon WW [27,43]. Our second code uses a different implementation of a multi- 
channel Monte Carlo generator with adaptive weight optimization. 

2.3 Classification of QCD corrections 

As QCD corrections to Higgs production via VBF we consider the interference of VBF 
diagrams of the type shown in Figure 1 with the virtual QCD corrections arising from 
gluon exchange, gluon fusion, and gluon splitting. We also take into account the contribu- 
tions from real gluon emission and gluon-induced processes. We classify these corrections 
in the same way as done for the QCD corrections to H — ^ 4/ described in Ref. [22] upon 
considering possible contributions to the squared amplitude. The ampUtude itself receives 
contributions from one of the two generic tree diagrams shown in Figure 2 or from both. 
Thus, the square of this amplitude receives contributions from squared and interference 
diagrams of the types depicted in Figure 3. Type (A) corresponds to the squares of each 
of the Born diagrams, type (B) to their interference if two Born diagrams exist. 
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Figure 4: Categories of interference diagrams contributing to the QCD corrections. 



After this prehminary consideration we define four different categories of QCD cor- 
rections. Examples of interference diagrams belonging to these categories are shown in 
Figure 4, the corresponding virtual QCD correction diagrams are depicted in Figure 5. 

(a) "Diagonal" QCD corrections to squared tree diagrams comprise all interference di- 
agrams resulting from diagram (A) of Figure 3 by adding one additional gluon. 
Cut diagrams in which the gluon does not cross the cut correspond to virtual one- 
loop corrections, the ones where the gluon crosses the cut correspond to real gluon 
radiation. Note that interference diagrams in which the gluon connects the two 
closed quark lines identically vanish, because their colour structure is proportional 
to Tr(A'')Tr(A'') = 0, where is a Gell-Mann matrix. Thus, the only relevant 
one-loop diagrams in this category are gluonic corrections to the Vqq' vertex, as 
illustrated in the first diagram of Figure 5; the real corrections are induced by the 
corresponding gluon bremsstrahlung diagrams. 

Previous calculations [11,16-19] of NLO QCD corrections focused on this category of 
corrections to t- and w-channel contributions only. This approximation is motivated 
by the smallness of s-channel contributions, at least in the kinematic domain relevant 
for Higgs production via VBF, and by the suppression of all types of interferences 
in lowest order. Both of these suppressions are due to strong enhancements in the 
t- and M-channel weak-boson propagators that receive a small momentum transfer; 
only in contributions to squared amplitudes that are related to squared t- and u- 
channel LO graphs four enhancement factors of this kind can accumulate. For 
instance, interferences between two different t- and w-channel tree diagrams involve 
four enhanced propagators, but they pairwise peak in different regions of phase 
space (forward or backward scattered quarks). 

(b) QCD corrections to interferences comprise all interference diagrams resulting from 
diagram (B) of Figure 3 by adding one additional gluon, analogously to the previous 
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(a,b) 



(b) 



H 




H 





(d) 



H 




H 




Figure 5: Basic diagrams contributing to the virtual QCD corrections to qq qqB. where 
V = W, Z and Q = d, u, s, c, b,t. The categories of QCD corrections, (a)-(d), to which 
the diagrams contribute are indicated. 

category. Relevant one-loop diagrams are, thus, vertex corrections or pentagon 
diagrams, as illustrated in the first two diagrams of Figure 5. 

(c) Corrections induced by one qqg splitting result from loop diagrams exemplified by 
the third graph in Figure 5. The remaining graphs are obtained by shifting the 
gluon to different positions at the same quark line and by interchanging the role of 
the two quark lines. Thus, the diagrams comprise not only box diagrams but also 
vertex diagrams. They do not interfere with Born diagrams with the same fermion- 
number flow because of the colour structure, i.e. in 0{as) they only contribute if 
two Born diagrams exist. 

Some of the squared diagrams of this category actually correspond to (collinear- 
singular) real NLO QCD corrections to loop-induced H-l-ljet production, e.g. qg — > 
gH or qq — »• g'R. Here we consider only the interference contributions of the loop 
diagrams of this category with the lowest-order diagrams where the Higgs boson 
couples to a weak boson (see Figure 1), resulting in a UV- and IR- (soft and collinear) 
finite correction. 

(d) Corrections induced by two qqg splittings (gg fusion) result from diagrams exempli- 
fied by the fourth graph in Figure 5. There are precisely two graphs with opposite 
fermion-number flow in the loop. Again, owing to the colour structure (see also 
below), these diagrams do not interfere with Born diagrams with the same fermion- 
number flow, i.e. the existence of two Born diagrams is needed. 

The squared diagrams of this category actually correspond to (coUinear-singular) 
real NNLO QCD corrections to loop-induced Higgs production via gluon fusion. 
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gg — s> H. The considered interference contributions of the loop diagrams of this 
category with the lowest-order diagrams of Figure 1, however, again yield a UV- 
and IR- (soft and coUinear) finite correction. 

This category of QCD corrections was recently considered in the approximation of 
an infinitely heavy top quark in Ref. [44] and found to be suppressed. There it was 
also argued that QCD corrections to these small contributions might be sizeable, 
because further gluon exchange between the two incoming (anti-) quarks enables 
an interference with the tree diagram with the same fermion-number fiow, thereby 
receiving an enhancement by four propagators with small momentum transfer. This 
contribution has very recently been studied in Ref. [45] and found to be completely 
negligible owing to the appearance of several other suppression mechanisms. 

2.4 Structure of virtual corrections 

Since the colour flow in EW loop diagrams is the same as in the corresponding lowest- 
order diagrams, the EW one-loop amplitudes ^ew can be decomposed into colour- and 
CKM-stripped amplitudes tMew exactly in the same way as done in lowest order, where 
we decomposed in terms of Ai^^ (2-2). 

According to the above classification, the QCD one-loop amplitudes of category (a) 
as well as the vertex corrections of category (b) involve only the colour operator Ci of 
(2.3), while the pentagon diagrams of category (b) and all loops of categories (c) and (d) 
involve only the colour operator C2. Thus, we can decompose the amplitudes ^QCD,(a)) 
etc., into colour- and CKM-stripped parts A^QCD(a)' ^^^-i ^ follows 

.Uidj^Ukdi _ T/; T/*^ciC244 . .WW (^oo'^'\ 

'^QCD(a)+QCD(fo,vert) ~ ^ij^kl'-'l •'^'QCD(a)+QCD(6,vert) 1"^^^ J 

- SikSjlCi^ ^ ^ ^•^QCD(a)+QCD(6,vert)(ll'2'2), 

Au^dj^Ukdi _ -r^ T/*^ciC244 . -WW f^oo'^'\ 

•^QCD(6,pent)+QCD(c)+QCD(ci) — ^ij^kl^2 •'^lQCD(&,pent) V-^^^ J 

- ^ik^jlC2^ ^ ^ ^•^QCD(6,pent)+QCD(c)+QCD(d)(ll'2'2), 

.didj^dkd, —A A r"=i'=2C2c; . ,zz 

•^QCD(a)+QCD(6,vert) ~ "ijOkl^l •'^'QCD(a)+QCD(6,vert) l-"-^^ -L ) 

- ^ik^jlCi ^ ^ ^-^QCD(a)+QCD(6,vert)(ll'2'2), 

.didj^dkdi —A A r"=i'^2C2c; . -zz (^oo'^'\ 

•^QCD(6,pent)+QCD(c)+QCD((i) " (^ij^kl^l •'^lQCD(6,pent)+QCD(c)+QCD(d) l"^^^ ^ ) 

- ^ik^jlC2^ ^ ^ ^•^QCD(6,pent)+QCD(c)H-QCD((i)(ll'2'2). (2.7) 

Note that W-mediated parts Al^^ do not receive contributions of categories (c) and (d). 

Since the lowest-order amplitudes only involve colour operators Ci, the following colour 
sums appear in the calculation of squared lowest-order amplitudes and of interferences 
between one-loop and lowest-order matrix elements: 

a,h,c,d a,h,c,d 

Xj^) ^ {Cf^*Cf'^) = N^, Xi^^ = {Cf"^*Cf'^) = N^, (2.8) 

a,b,c,d a,b,c,d 
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where Y,a,b,c,d stands for the sum over the colour indices a, b, c, d, and iV^ = 3 is the colour 
factor for a quark. Squared Born diagrams, as illustrated in type (A) of Figure 2, are 
proportional to x['^\ lowest-order interference diagrams of type (B) are proportional to 
x[^\ The situation is analogous for all EW one- loop diagrams. By definition, category (a) 
of the gluonic diagrams comprises all one-loop QCD corrections proportional to X^'^K In 
category (b), the vertex corrections are proportional to x[^^ and the pentagons to X2^\ 
Categories (c) and (d) receive only contributions from X2^^; interferences of one-loop 
diagrams like (c) and (d) in Figure 5 with Born diagrams of the same fermion-numbcr 
flow vanish because of X^^^ = 0. Finally, for the one-loop corrections to the squared 
matrix elements we obtain 

+ 5.fc5,,>io^^(ll'2'2)*[9A^i^+QCD(a)(ll'2'2) - 3l-.,l-felA^^^+QCD(;.)(122'l')] }, 



^didj^dkdi 



J2 2Re|(^f''^^^"'=*)*^- 
{ci,<} 

= 2Re{5,,-4iMf(122'l')*[9M|VQCD(a)(122'l') 

~35ifc5j7A^|w+QCD(6)+QCD(c)+QCD{d)(ll'2'2) 



+ 5,fc5,z-Mo''(ll'2'2)* [9>t^^+QCD(a)(ll'2'2) 

-3%4/-MEW+QCD{6)+QCD(c)+QCD(d)(122'l')j |. (2.9) 

As already observed for the squared LO amplitudes, also here only absolute values of the 
CKM matrix elements, such as jV^jp, contribute after contracting the Kronecker deltas 
of the generation indices. 

2.5 Hadronic cross section 

The hadronic cross section (Tpp(Pi,P2) for colliding protons results from the partonic 
cross section 0"jj(xiPi, ^2^2) upon convolution with the parton densities fi{xi,fip), which 
corresponds to parton i carrying the fraction xi of the proton momentum P; (/ = 1, 2), 

(7pp(Pi,P2) = / dXi [ dX2^fi{xi,f4)fj{x2,t4)aij{xiPi,X2P2), (2.10) 

Jo Jo 

where /xp is the factorization scale that separates the hard partonic process from the 
soft physics contained in the PDFs. The sum over the partons i,j includes all quarks, 
antiquarks, gluons, and the photon. In LO only quarks and/or antiquarks are present in 
the initial state, in NLO also processes with one gluon or photon contribute. In detail the 
NLO parton cross sections read 
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(J, 



97 



'^97, real ~l~ '^g7,fact; ''"7'? '-'"79, real ^79, fact) 



cr, 



g7 



cr^ 



cr. 
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0, 



(2.11) 



where q generically stands for any relevant quark or antiquark. The LO and virtual one- 
loop contributions ("LO" and "virt") involve the partonic 2^2 kinematics, while real 
emission contributions ("real") are of the type 2 — > 3 with one additional light (anti)quark, 
gluon, or photon in the final state. The calculation of these subcontributions has been 
briefly described in the previous sections. The contribution called "fact" results from the 
PDF redefinition necessary to absorb collinear initial-state singularities into the PDFs via 
factorization, so that the partonic cross sections (Xjj are free of such singularities. This 
separation introduces a logarithmic dependence on the factorization scale /ip in ajj^fact that 
compensates the implicit //p dependence in the PDFs in NLO accuracy. The factorization 
explicitly proceeds as follows. 

The virtual and real contributions of the parton cross sections contain mass singulari- 
ties of the form «s In(mg) and a In(mq), which are due to collinear gluon/photon radiation 
off the initial-state quarks or due to a collinear splitting g/7 qq of initial-state glu- 
ons or photons. For processes that in LO involve only quarks and/or antiquarks in the 
initial state, the factorization is achieved by replacing the (anti-)quark distribution fq{x) 
according to (see e.g. Ref. [41]) 



fq{x,l^l) ^ fq{x,lll) - f —fq(-,l4 

Jx Z \Z 



)(21n(l-2) + l)] +C 



//,QCDl 



^^//(^)(21n(l-^) + l)]_^ + C// ,qed(^) I f 



(2.12) 



where Cij{z) are the so-called coefficient functions, and the splitting functions Pij{z) are 
defined as 



1 + z^ 
1-z 



Pf,{z) = z' + {l 



(2.13) 



Starting from the hadronic LO cross section after the substitution (2.12), the factorization 
contributions o"jj,fact correspond to the terms of 0{as) and 0{a) involving the PDF of i 
and j. The replacement (2.12) defines the same finite coefficient functions as the usual 
D-dimensional regularization for exactly massless partons where the In(mq) terms appear 
as 1/(D — 4) poles. The actual form of the coefficient functions defines the finite parts of 
the NLO corrections and, thus, the factorization scheme. Following standard definitions 
of QCD, we distinguish the MS and DIS-like schemes which are formally defined by 



Cf^iz) = Cl'{z) = C)f (.) = 0, 



/7 
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C]?f(.) 



47 4 



C^f (z) = C^f (z) = Pj^(z) In(^) - ^z" + 8z - 1. (2.14) 

The MS scheme is motivated by formal simphcity, because it merely rearranges the IR- 
divergent terms (plus some trivial constants) as defined in dimensional regularization. The 
DIS-like scheme is defined in such a way that the deep inelastic scattering (DIS) structure 
function F2 does not receive any corrections; in other words, the radiative corrections to 
electron-proton DIS are implicitly contained in the PDFs. 

Whatever scheme has been adopted in the extraction of PDFs from experimental data, 
the same scheme has to be used when predictions for other experiments are made using 
these PDFs. In particular, the absorption of the collincar singularities of both QCD 
and QED origin into PDFs requires the inclusion of the corresponding QCD and QED 
corrections into the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution of 
these distributions and into their fit to experimental data. We use the MRST2004QED 
PDFs [46] which consistently include QCD and QED NLO corrections. These PDFs 
include a photon distribution function for the proton and thus allow to take into account 
photon-induced partonic processes. As explained in Ref. [41], the consistent use of these 
PDFs requires the MS factorization scheme for the QCD corrections, but the DIS scheme 
for the QED corrections, i.e. we employ Cff and of (2.14) for the QCD, but C^f 
and C^^^ for the QED corrections. 

3 Numerical results 

3.1 Input parameters and setup 

We use the following set of input parameters [47] , 

0.1187, 



— 


1.16637 X 10-5 GeV-^ 


a{0) = 


1/137.03599911, 


as(Mz) 


/M'LEP_ 


80.425 GeV, 


pLEP 
J- w ~ 


2.124 GeV, 




M^EP = 


91.1876 GeV, 


pLEP 

J- z — 


2.4952 GeV, 




mc = 


0.51099892 MeV, 


mf,= 


105.658369 MeV, 


rrir 


m^ = 


66 MeV, 


mc = 


1.2 GeV, 


mt 


md = 


66 MeV, 


ms = 


150 MeV, 





1.77699 GeV, 
174.3 GeV, 
4.3 GeV. (3.1) 

If not stated otherwise, the Higgs-boson mass is set to 

MH = 120GeV. (3.2) 

Using the complex- mass scheme [34], we employ a fixed width in the resonant W- and Z- 
boson propagators in contrast to the approach used at LEP to fit the W and Z resonances, 
where running widths are taken. Therefore, we have to convert the "on-shell" values of 
M^^P and T\.^^ {V = W,Z), resulting from LEP, to the "pole values" denoted by My 
and Fy. The relation between the two sets of values is given by [48] 



Mv = M^EP/ Ji + (rLEP/MLEP)2^ Yy = T\P'' / J 1 + {T\.^^ / Ml^^^ ^ , (3.3) 
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leading to 



Mw = 80.397 ... GeV, Tw^ 2.123 ... GeV, 
Mz= 91.1535... GeV, Tz = 2.4943 ... GeV. (3.4) 

We make use of these mass parameters in the numerics discussed below, although the 
difference between using My or My^^ would be hardly visible. 

The masses of the light quarks are adjusted to reproduce the hadronic contribution to 
the photonic vacuum polarization of Ref. [49]. Since quark mixing effects are suppressed^ 
we neglect quark mixing and use a unit CKM matrix. 

We use the scheme, i.e. we derive the electromagnetic couphng constant from the 
Fermi constant according to 

ac, = V2G^MUi - M^/MD/tt. (3.5) 

In this scheme, the weak corrections to muon decay Ar are included in the charge renor- 
malization constant (see e.g. Ref. [50]). As a consequence, the EW corrections are prac- 
tically independent of the masses of the hght quarks. Moreover, this definition effectively 
resums the contributions associated with the running of a from zero to the W-boson mass 
and absorbs leading universal corrections oc G^m\ from the p parameter into the LO 
amplitude. 

We use the MRST2004QED PDFs [46] which consistently include 0{a) QED correc- 
tions. Since no associated LO PDFs exist, we use these distributions both for LO and 
NLO predictions. We do not include processes with external bottom quarks in our default 
set-up. These are suppressed either because of the smallness of the b-quark densities or 
due to s-channel suppression. Partonic processes involving b quarks are, however, in- 
cluded in our code in LO. As discussed in Section 3.4, these contributions are at the level 
of a few per cent. In contrast to Ref. [20], we use Mw (instead of Mh) as factorization 
scale both for QCD and QED coUinear contributions, which is a better scale choice when 
considering large Higgs-boson masses. For the calculation of the strong coupling constant 
we employ Mw as the default renormalization scale, include 5 flavours in the two-loop 
running, and flx a^^Mz) = 0.1187. 

Jet reconstruction from final-state partons is performed using the /cx-algorithm [51] as 
described in Ref. [52]. Jets are reconstructed from partons of pseudorapidity l^yj < 5 using 
a jet resolution parameter D — 0.8. Real photons are recombined with jets according to 
the same algorithm. Thus, in real photon radiation events, final states may consist of jets 
plus a real identifiable photon, or of jets only. 

We study total cross sections and cross sections for the set of experimental "VBF 
cuts" defined in Ref. [18]. These cuts are expected to significantly suppress backgrounds 
to VBF processes, enhancing the signal-to-background ratio. We require at least two hard 
jets with 

PTj>20GeV, bj|<4.5, (3.6) 

where pxj is the transverse momentum of the jet and yj its rapidity. Two tagging jets ji and 
j2 are defined as the two jets passing the cuts (3.6) with highest pt such that pxji > PTj2- 

^We checked that the cross section without cuts changes by one per mille and the one with VBF cuts 
by less than 0.01% when using a realistic quark mixing matrix. 
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Furthermore, we require that the tagging jets have a large rapidity separation and reside 
in opposite detector hemispheres: 

^Vii = IVh - %2l > 4, Z/ji • < 0. (3.7) 

All presented results have been obtained using the subtraction method. For the results 
in the tables we used 10^ events for the setup with VBF cuts and 5 x 10^ events without 
cuts. For the plots of Mh and factorization scale dependence we generated 10^ events 
without cuts and 2 x 10^ events with VBF cuts. The plots for the distributions are based 
on 10^ events. Generally, the real corrections and the finite virtual QCD corrections are 
only calculated for each 10th event, the finite virtual EW corrections only for each 100th 
event. 

3.2 Results for integrated cross sections 

We first consider results for integrated cross sections. In Figure 6 we plot the total 
cross section with and without VBF cuts as a function of the Higgs-boson mass. In the 
left panel we show the absolute predictions in LO and in NLO including QCD and EW 
corrections. For Mh = 100 GeV to about 200 GeV the results without cuts are larger by a 
factor 2-4, while for Mh = 700 GeV this factor reduces to 1.7. In the right panel we show 
the relative QCD and EW corrections separately. For Higgs-boson masses in the range 
100-200 GeV, without cuts, the QCD corrections drop from +5% to 0%, and the EW 
corrections are about —5% depending only weakly on the Higgs-boson mass. For very 
small Higgs-boson masses, QCD and EW corrections cancel each other substantially. With 
VBF cuts the EW corrections are somewhat more negative, while the QCD corrections 
vary between —4% and —6%. For higher Higgs-boson masses, the QCD corrections do 
not change much and reach 1% and —7% at Mh = 700 GeV without cuts and with VBF 
cuts, respectively. The EW corrections increase steadily with the Higgs-boson mass up to 
8% and 7% at Mh = 700 GcV without cuts and with VBF cuts, respectively. In the EW 
corrections the WW, ZZ, and tt thresholds are clearly visible. 

It is interesting to note that, at least for Higgs-boson masses below 200 GeV, the 
EW corrections to the full VBF channel are similar in size and sign to the subreactions 
pp ^ WH/ZH + X [50]. Compared to the related decays H ^ WW/ZZ ^ 4/ [21,22] the 
size is similar, but for low Higgs masses of 100— 200 GeV the sign is different. 

In Table 1 we present numbers for integrated cross sections for Mh = 120, 150, 200, 
400, and 700 GeV without any cuts and in Table 2 results for the VBF cuts defined 
above. We hst the LO cross section ctlo, the cross section including NLO QCD and EW 
corrections, (Jnlo, and various contributions to the relative corrections. The complete 
EW corrections Sew comprise the EW corrections resulting from loop diagrams and real 
photon radiation, SEw,qq, and the corrections from photon-induced processes Sq.y. Fur- 
thermore, 5-Ew,qq includes the dominant two-loop correction 5q2 due to Higgs-boson 
self- interaction, which was introduced in Section 2.2. The QCD corrections 5qcd are de- 
composed in the diagonal contributions ^qcD.diag, non-diagonal contributions (^QCD,nondiag, 
the contribution resulting from gluon splitting 5g-spiit, and those from gluon-gluon fusion 
(5gg_fusion, as explained in Section 2.3. 

The QCD corrections are dominated by the diagonal contributions, i.e. by the vector- 
boson-quark-antiquark vertex corrections to squared LO diagrams. All other contribu- 
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Figure 6: Higgs-boson-mass dependence of LO and complete NLO cross section (left) and 
relative EW and QCD corrections (right) without and with VBF cuts. 

tions are at the per-mille level and even partially cancel each other. They are not enhanced 
by contributions of two t- or -u-channel vector bosons with small virtuality and therefore 
even further suppressed when applying VBF cuts. The photon-induced EW corrections 
are about 1% and reduce the EW corrections for small and intermediate Mh. The two- 
loop correction ^g^m* is negligible in the 1ow-Mh region, but becomes important for large 
Higgs-boson masses. For Mh = 700 GeV this contribution yields +4% and constitutes 
about 50% of the total EW corrections. Obviously for Higgs masses in this region and 
above the perturbative expansion breaks down, and the two-loop factor 5q2 might serve 
as an estimate of the theoretical uncertainty. 

3.3 Subcontributions from s channel and t/u interference 

Previous calculations of the VBF process [11,16-19] have consistently neglected s- 
channel contributions ( "Higgs strahlung" ) , which involve diagrams where one of the vector 
bosons can become resonant, as well as the interference between t- and w-channel fusion 
diagrams. To better understand the effect of these approximations we have calculated 
these contributions to the integrated cross section. In Table 3 and Table 4 we present, 
with and without VBF cuts, respectively, contributions from s-channel processes, cXs, and 
from t/w-channel interference terms (Tt/u-mt) at both LO and NLO. The NLO result does 
not include the corrections due to photon-induced processes, which cannot be split into 
the above subcontributions respecting gauge invariance. 

While t /^-interference terms, with or without VBF cuts, contribute less than 1% to 
the cross section, s-channel contributions are clearly non-negligible when no cuts are used. 
At LO (NLO), for Mh = 120 GeV they contribute 22% (27%) to the total cross section, 
while for Mh = 200 GeV this contribution decreases to 9% (10%). For Mh = 700 GeV 
s-channel processes contribute less than 1% to the cross section, with and without VBF 
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'^gg-fusion [%] 


-0.19(1) 


-0.20 


-0.22 


-0.24 


-0.11(1) 


5g2M4 [%] 


0.0027 


0.0073 


0.025 


0.42 


4.03(1) 



Table 1: Cross section for pp — > H + 2jets + X in LO and NLO without cuts and relative 
EW and QCD corrections split into various subcontributions. 



Mh [GeV] 


120 


150 


200 


400 


700 


o-LO [fb] 


1876.3(5) 


1589.8(4) 


1221.1(3) 


487.31(9) 


160.67(2) 


O-NLO [fb] 


1665(1) 


1407.5(8) 


1091.3(5) 


435.4(2) 


160.36(5) 


Sew [%] 


-6.47(2) 


-6.27(2) 


-4.98(1) 


-3.99(1) 


6.99(2) 


SEW,qq [%] 


-7.57(2) 


-7.42(2) 


-6.19(1) 


-5.37(1) 


5.44(2) 


Sq, [%] 


1.10 


1.15 


1.22 


1.38 


1.55 


Sqcd [%] 


-4.77(4) 


-5.20(4) 


-5.65(3) 


-6.67(3) 


-7.18(2) 


f^QCD.diag [%] 


-4.75(4) 


-5.17(4) 


-5.66(4) 


-6.63(3) 


-7.18(2) 


(^QCD.nondiag [%] 


-0.011 


-0.0052(1) 


0.0032(1) 


0.0030 


0.0022 


(^g-split [%] 


-0.0085(1) 


0.0084(1) 


0.027 


0.014 


0.0074 


'^gg-fusion [%] 


-0.030 


-0.030 


-0.028(1) 


-0.020 


-0.014 


Sgim^ [%] 


0.0035 


0.0086(1) 


0.027 


0.43 


4.06(1) 



Table 2: As in Table 1, but with VBF cuts applied. 



cuts. Thus, for increasing Higgs-boson masses, the contribution from Higgs-strahlung 
processes becomes less and less important compared to the contribution from pure fusion 
processes. When VBF cuts are used, both the s-channel and t /w-channel-interference 
contributions are strongly suppressed, yielding less than 0.6% of the cross section for all 
the studied Higgs-boson masses. The comparably large NLO s-channel contribution after 
VBF cuts originates from real gluon corrections with up to three jets in the final state, 
because in contrast to LO the two jets from the weak-boson decay, which tend to be 
aligned owing to a boost, are not forced to be the two well-separated tagging jets. We 
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0'LO,t/«-int [fb] 


-9.2 


-5.6 


-2.71 


-0.32 


-0.041 


CrNLO,t/«-int [fb] 


-27.6 


-9.4 


0.04(1) 


-1.08(1) 


-0.19 



Table 3: s-channcl contributions and contributions from interference between t and u 
channels to the pp — > H + 2jets + X cross section at LO and NLO, without any cuts. 



Mh [GeV] 


120 


150 


200 


400 


700 


o-LO,s [fb] 


0.0025 


0.0015 


0.00071 


0.000072 


0.0000069 


C^NLO.s [fb] 


9.45(1) 


5.21(1) 


2.33 


0.29 


0.044 


C^LO.t/u-int [fb] 


-0.12 


-0.091 


-0.060 


-0.016 


-0.0034 


<7NLO,t/u-int [fb] 


-0.75 


0.17 


0.76 


0.089 


0.0044(1) 



Table 4: As in Table 3, but with VBF cuts applied. 



conclude that, applying typical experimental VBF cuts, the contributions from s-channel 
diagrams and i/u-channel interferences can be safely neglected. 

3.4 Leading-order b-quark contributions 

In this section we present the contributions arising at LO from processes that include 
b-quarks in the initial and/or final states. There are three types of contributions involving 
b quarks. The first type consists of s-channel diagrams with a bb pair in the initial state, 
the second type comprises s-channel diagrams with a bb pair in the final state, and the 
third type involves s-channel diagrams with a bb pair in both the initial and the final 
state as well as all t- and w-channel diagrams where a b or b quark goes from the initial 
state to the final state. In Table 5 we show, for different Mh values, LO cross-section 
results without b-quark contributions, (Jlo, nob, the results including only initial-state 
b quarks, ctlo, b-in, the results including only final-state b quarks, ctlo, b-out, and including 
both initial- and final-state b quarks, ctlo, b-in/out- The relative contributions arising from 
these subprocesses, 5b-in, ^b-out, and 5b-in/out; are also shown. In Table 6 we present LO 
results including VBF cuts. 

For low Higgs-boson masses and no cuts, final-state b quarks increase the LO cross 
section by up to 2%. The increase due to initial-state b quarks is one per mille or less, being 
strongly suppressed due to the two bottom densities involved (bb-annihilation processes). 
Including both initial and final-state b quarks increases the total cross section by up to 4%, 
a contribution that is similar in absolute value to the total EW correction, but opposite 
in sign. The contributions from final-state and/or initial-state b quarks decrease with 
increasing Higgs-boson mass. For Mh = 700 GeV, b-quark contributions from cither final 
state or initial state become negligible, while simultaneous initial- and final-state b-quark 
corrections decrease to 1.4%. When VBF cuts are imposed, b-quark contributions become 
less important. This is particularly noticeable in contributions arising from processes with 
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0"LO, b-in/out [ fb] 


6203(1) 


4495(1) 


2945.7(6) 


919.5(2) 


274.49(4) 


(^b-in/out [%] 


4.37(2) 


3.79(2) 


3.17(2) 


2.09(2) 


1.44(2) 



Table 5: LO cross section for pp — > H + 2jets + X with and without initial- and/or 
final-state b quarks, without any cuts. 



Mh [GeV] 


120 


150 


200 


400 


700 


O-LO, nob [fb] 


1876.1(5) 


1589.8(4) 


1221.1(3) 


487.32(9) 


160.66(2) 


0-LO, b-in [ fb] 


1876.1(5) 


1589.8(4) 


1221.1(3) 


487.32(9) 


160.66(2) 


5b-i„ [%] 


0.00(3) 


0.00(2) 


0.00(2) 


0.00(2) 


0.00(1) 


0-LO, b-out [fb] 


1876.1(5) 


1589.8(4) 


1221.1(3) 


487.32(9) 


160.66(2) 


<5b-out [%] 


0.00(3) 


0.00(2) 


0.00(2) 


0.00(2) 


0.00(1) 


0"LO, b-in/out [fb] 


1918.5(5) 


1624.5(4) 


1246.3(3) 


495.55(9) 


162.75(2) 


^b-in/out [%] 


2.26(3) 


2.18(2) 


2.06(2) 


1.69(2) 


1.30(1) 



Table 6: As in Table 5, but with VBF cuts applied. 



final-state but no initial-state b quarks and vice versa, Higgs-strahlung processes of the 
form qq — > bbif and bb — > qqH. These are s-channel processes and, as already shown in 
Section 3.3, this type of contributions are strongly suppressed by the VBF cuts. 

3.5 Scale dependence 

In Figures 7 and 8 we show the dependence of the total cross section on the factoriza- 
tion and renormalization scale for Mh = 200 GeV and Mh = 400 GeV, respectively. We 
relate the factorization scale /ip, which applies to both QCD and QED contributions, and 
the renormalization scale /xr to the W-boson mass as 

/X = //F = ^fMw, A«r = CrMw, (3.8) 

and vary and between 1/8 and 8. We study the scale dependence of the LO 
cross section, of the QCD-corrected NLO cross section, and of the complete NLO cross 
section including both QCD and EW corrections for ^r = ^p- In addition wc depict 
the QCD-correctcd NLO cross section for the setup where ^r = 1/^f (NLO QCD'). For 
Mh = 200 GeV, varying the scale up and down by a factor of 2 (8) changes the cross 
section by ±3% (±9%) in LO and by ±1% (-2%/+9%) in NLO for the set-up without 
cuts. With VBF cuts, the scale uncertainty amounts to ±6% (±18%) in LO and ±2% 
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Figure 7: Scale dependence of LO and NLO cross section with QCD or QCD+EW cor- 
rections for Mh = 200 GeV without cuts (left) and with VBF cuts (right). 
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Figure 8: Scale dependence of LO and NLO cross section with QCD or QCD+EW cor- 
rections for Mh = 400 GeV without cuts (left) and with VBF cuts (right). 
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(±11%) in NLO. For Mr = 400 GeV, the scale uncertainty is reduced from ±7% (±24%) 
in LO to ±1% (±11%) in NLO for the cross section without cuts, and from ±8% (±29%) 
in LO to ±3% (±15%) in NLO for the cross section with VBF cuts. For Mh = 400 GeV, 
it is clearly seen from the results that Mw is a more appropriate scale choice than Mh. 
For this reason we have chosen Mw as default scale in this paper, while we used Mh in 
Ref. [20], where we only considered Higgs-boson masses comparable to Mw- 

3.6 Slicing cut dependence 

In the slicing approach (as e.g. reviewed in Ref. [53]), phase-space regions where real 
photon/gluon emission and photon/gluon-induced processes contain soft or coUinear sin- 
gularities are defined by the auxiliary cutoff parameters Ss, 5c <^ 1- In real photon/gluon 
radiation processes, the region 

X<k'< S,^, (3.9) 

where k is the photon/gluon momentum, the partonic centre-of-mass energy, and A an 
infinitesimal photon/gluon mass, is treated in soft approximation. The regions determined 

by 

1 - cos(%g},) < 5e, > S,—, (3.10) 

where O^^^g^g is the angle between any quark q and the photon or gluon, are evaluated using 
coUinear factorization. In photon- or gluon-induced processes, singularities arise only in 
the coUinear region, i.e. a slicing cut on the angle between any final-state quark q and 
the initial-state photon or gluon is sufficient to exclude the singularity from phase space. 
Specifically, we define this angular cut as in (3.10). The coUinear-splitting singularities 
are also treated using coUinear factorization. 

In the remaining phase space no regulators (photon/gluon and quark masses) are used. 
Therefore, the slicing result is correct up to terms of 0{Ss) and 0{Sc). In Figure 9 we 
show the dependence of the complete corrections to the cross section with VBF cuts on 
Ss for fixed 6c = 10~^ and the dependence on 6c for fixed 6s = 10^^. The error bars rcficct 
the uncertainty of the Monte Carlo integration. These results were obtained with lO'^ 
events for the slicing method and 10^ events for the subtraction method, using Mh as 
factorization and renormalization scale. For decreasing auxiliary parameters Ss and 5c, 
the slicing result reaches a plateau and becomes compatible with the subtraction result. 
The integration error in the result obtained with the slicing method increases for lower 
cut-off parameters. On the other hand, the subtraction results, for the same number of 
events, always show smaller integration errors. 

3.7 Differential cross sections 

In this section we consider results for distributions involving Higgs-boson and tagging- 
jet observables. We show results for Mh = 120 GeV in the setup including VBF cuts. For 
each distribution we plot the absolute predictions in LO and in NLO including QCD and 
EW corrections. In addition, we show the relative corrections, both the QCD and EW 
corrections separately, as well as their sum. 
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Figure 9: Dependence of the relative corrections to the total cross section with VBF cuts 
on the energy cutoff 4 for 5c = 10~^ (l.h.s.) and on the angular cutoff 6c for 6s = 10~^ 
(r.h.s.) in the slicing approach for Mh = 120 GeV. For comparison the corresponding 
result obtained with the dipole subtraction method (with 10 times less statistics) is shown 
as a la band in the plots. 



We first consider Higgs-boson observables and show the distribution in the transverse 
momentum pt,h in Figure 10. The differential cross section drops strongly with increasing 
Pt,h, while both the relative EW and QCD corrections increase in size and reach —20% 
for pt,h = 500 GeV. It is interesting to note the differences between this result and the 
same distribution in Higgs-boson production via gluon-fusion, as e.g. shown in Ref. [54]. 
In weak-boson fusion, this distribution is broader and peaks at a much larger value of 
Pt,h- 

The distribution in the rapidity ?/h of the Higgs boson is presented in Figure 11. While 
the relative EW corrections depend only weakly on this variable, the QCD corrections 
show an increase for large rapidities. Total corrections decrease the differential cross 
section by more than 10% in the central region, inducing an important change in the 
shape of this distribution. 

Figures 12 and 13 show the differential cross section as function of the transverse 
momentum of the harder and softer tagging jet, respectively. These distributions peak 
near or below pj^x ~ and then drop strongly with increasing jet transverse momentum. 
QCD and EW corrections become more and more negative with increasing pj^T- For low 
transverse momentum these corrections are at the level of 5%, while for p^^x = 400 GeV 
and 150 GeV they add up to about —38% and —25% for the harder and softer tagging 
jet, respectively. This induces a substantial change in shape of these distributions. 

In Figures 14 and 15, we depict the distributions in the rapidities of the harder and 
softer tagging jet, respectively. It can be clearly seen that the tagging jets are forward 
and backward located. The EW corrections vary between —4% and —7%. The QCD 
corrections exhibit a strong dependence on the jet rapidities. For the harder tagging jet 
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Figure 10: Distribution in the transverse momentum pt,h of the Higgs boson (left) and 
corresponding relative corrections (right) for Mh = 120 GeV. 
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Figure 11: Distribution in the rapidity ?/h of the Higgs boson (left) and corresponding 
relative corrections (right) for Mh = 120 GeV. 
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Figure 12: Distribution in the transverse momentum Pji,T of the harder tagging jet (left) 
and corresponding relative corrections (right) for Mh = 120 GeV. 
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Figure 13: Distribution in the transverse momentum Pj2,T of the softer tagging jet (left) 
and corresponding relative corrections (right) for Mh = 120 GeV. 
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Figure 14: Distribution in the rapidity y^^ of the harder tagging jet (left) and correspond- 
ing relative corrections (right) for Mh = 120 GeV. 
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Figure 15: Distribution in the rapidity y^^ of the softer tagging jet (left) and corresponding 
relative corrections (right) for Mh = 120 GeV. 

they are about —8% in the central region but become positive for large rapidities, where 
they tend to compensate the EW corrections. For the softer tagging jet the variation for 
large rapidities is smaller, and the QCD corrections become small also near yj = 0. Shape 
changes due to the full corrections can reach 10%. 

In Figure 16 we present the distribution in the azimuthal angle separation of the two 
tagging jets. This distribution is particularly sensitive to non-standard contributions to 
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Figure 16: Distribution in the azimuthal angle difference A0jj of tlie tagging jets (left) 
and corresponding relative corrections (right) for Mh = 120 GeV. 
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Figure 17: Distribution in the tagging-jet-pair invariant mass Mjj (left) and corresponding 
relative corrections (right) for Mh = 120 GeV. 

the vertices [18]. As expected for VBF processes, there is a large azimuthal angle 

separation between the two tagging jets. While QCD corrections are almost fiat in this 
variable, the QCD+EW corrections exhibit a dependence on A^jj on the level of 4%. 

Finally, in Figure 17 we show the distribution in the tagging-jet-pair invariant mass 
Mjj. Tagging jets identified in EW processes have typically larger jet-pair invariant masses 
than the ones identified in QCD processes. Consequently, Mjj can be used to further 
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SI 7 f^(^ ^ 
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^NLO [fb] 

<LO [fb] 


4494 f 41 
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3400(3) 


^oyjo 

2503.4(6) 
2666(2) 


R^R A(7\ 

OtJO.^l 1 I 

858.8(2) 
839.0(7) 


268.03(6) 
285.9(3) 



Table 7: Total cross section for pp — > H + 2jets + X in LO and NLO calculated with our 
program, ctlo/nlO; and with VV2H, Clq/nlO' f'^^' setup defined in the text. 



suppress QCD backgrounds, as e.g. pointed out in Ref. [1]. This distribution peaks at 
approximately Mjj = 500 GeV and is strongly suppressed for higher invariant-mass values. 
The EW corrections decrease with increasing Mjj and compensate the increasing QCD 
corrections for large invariant masses. The total correction is of the order of —10%. 

3.8 Comparison with related NLO QCD calculations 

In this section we compare our results to those obtained with the software packages 
VV2H by M. Spira [55] and VBFNLO by D. Zeppenfeld et al.[56]^. These programs allow 
to calculate the LO and NLO-QCD-corrected cross sections for Higgs-boson production 
via VBF at hadron colliders. It is important to note that s-channel contributions and 
i/w-channel- interference contributions are not taken into account in these calculations. In 
particular, only the 0{as) corrections that correspond to our class (a) of QCD contribu- 
tions (see Section 2.3) are included. In order to allow for a tuned comparison, we here use 
only four quark flavours for the external partons, i.e. we have switched off the effect of 
initial- and final-state b quarks in the calculations. We compare the results of VV2H and 
VBFNLO with LO and NLO-QCD-corrcctcd results of our code with s-channcl contribu- 
tions and t/-u-channel-interference contributions switched off, o'*"°°'^. In addition, we give 
the results of our code, a^"^^^, with these contributions and all interferences switched on 
and including all EW corrections apart from photon-induced processes. We use CTEQ6 
parton distributions [57] and our default set of input parameters. 

We first compare our results to those obtained with VV2H, which implements the 
formulae presented in Ref. [16]. As it is not possible to include phase-space cuts in VV2H, 
we have only compared total cross sections. The results of this comparison can be found 
in Table 7. We observe that the LO cross sections agree within 0.05% and the NLO 
corrected results within 0.2%, a difference which is of the order of the statistical error. 
Our complete predictions a^'^^^ differ from the results of VV2H by up to 30% for low 
Higgs-boson masses and by a few per cent for high Higgs-boson masses. The bulk of this 
big difference for small Mh values is due to the missing s-channel contributions in VV2H. 

We now turn to VBFNLO, which implements the results of Ref. [17]. As explained 
there, VBFNLO generates an isotropic Higgs-boson decay into two massless "leptons" (which 
represent t"^t~ or 77 or bb final states), and imposes a cut on the invariant mass of the 

^For this comparison we have employed the VV2H version dated July 23 2007, and VBFNLO-v .1.0. 
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Higgs boson. In order to be able to compare with this setup, we have implemented a 
convolution with a Breit-Wigner distribution for the Higgs-boson in one of our codes. 
When performing this convolution we can either evaluate the matrix element for Higgs 
production for an on-shell Higgs boson or for a Higgs boson with an invariant mass given by 
the Breit-Wigner distribution. While the first variant is gauge invariant and corresponds 
to a pole approximation, the second one, which is implemented in VBFNLO, violates EW 
gauge invariance. Because of the simple structure of the matrix element, this might not 
be a problem in LO and if only QCD corrections are included. Both variants neglect 
contributions that do not involve a resonant Higgs boson, which is a good approximation 
for small Higgs-boson masses, where the Higgs-boson width is small, but not for large 
Higgs-boson masses, where the width is large. Using these two variants of our code, we 
have compared cross sections without imposing any cuts on the decay products of the 
Higgs boson and using a unit branching ratio. To define the integration region in the 
neighbourhood of the Higgs resonance, we employ the value for the Higgs-boson width 
calculated by VBFNLO. 

The results for the cross section without cuts are compared in Table 8, while results 
including VBF cuts can be found in Table 9. The relative difference between the results 
of VBFNLO and the variant of our code with off-shell matrix elements, u*"'^^'^, is below 
0.04% for the total LO cross section and below 0.2% for the NLO-QCD-corrected cross 
section, both with and without VBF cuts. This difference is of the order of the statistical 
error. The difference between 0"*"'^'='^ and the variant with on-shell matrix elements, 
is at the per-mille level for Higgs-boson masses below 200 GeV but strongly increases for a 
heavy Higgs boson. For Mh = 400 GeV and 700 GeV the differences reach about 4% and 
30%, respectively, which illustrates the order of uncertainty without a more sophisticated 
treatment of off-shell effects of the Higgs boson including its decay. The results for u^*''^* 
are obtained with on-shell matrix elements only, since the off-shell matrix elements with 
EW corrections become gauge dependent. For small Higgs-boson masses and VBF cuts 
applied these predictions differ from those of VBFNLO by one per mille or less in LO and 
by 6-8%, the size of the EW corrections, in NLO. On the other hand, without cuts the 
big difference between a^*^^^ and the other predictions at small Mh values is again due to 
s-channel contributions. 

4 Conclusions 

Higgs-boson production via weak-boson fusion is one of the most important processes 
in the search for and the study of a Standard Model-like Higgs boson at the LHC. In this 
paper we present the first calculation of the NLO electroweak corrections for this process 
and we extend previously existing approximate NLO QCD calculations by including s- 
channcl topologies (Higgs-strahlung processes) and all interferences, both in LO and NLO. 

We find that the electroweak corrections are of the order of 5-10%, i.e. as large as the 
NLO QCD corrections. Real corrections induced by photons in the initial state increase 
LO results by roughly 1%. More precisely, the electroweak corrections are approximately 
—5% for Higgs masses below 200 GeV and for larger Mh values steadily increase up to 
about -1-7% for Mh = 700 GeV. For this Higgs-boson mass the leading two-loop effects 
in the heavy- Higgs limit, which are included in our calculation, become as large as the 
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Table 8: Cross section for pp ^ H + 2jets + X in LO and NLO calculated with our 
program, ctlo/nlo, and with VBFNLO, ctlq/nlO' without any cuts and for the setup defined 
in the text. 
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Table 9: As in Table 8, but with VBF cuts applied. 



one-loop corrections. This signals the breakdown of perturbation theory for large Higgs- 
boson masses. We suggest that the theoretical uncertainty from missing higher-order 
corrections can be estimated by the size of the leading two-loop heavy-Higgs effects in 
this domain. Moreover, for Mh ^ 400 GeV owing to the large Higgs-boson width the 
on-shell approximation is not sufficient any more and a more sophisticated treatment 
including off-shell effects of the Higgs boson and its decay width is required. 

We have implemented our calculation in a flexible Monte Carlo event generator, and 
studied differential distribution in Higgs-boson and tagging-jet observables. Specifically, 
we have presented results for distributions in transverse momenta, in rapidities, in the 
azimuthal angle difference of the tagging jets, and in the tagging-jet pair invariant mass. 
We found that QCD and electroweak corrections do not simply rescale differential distri- 
butions, but induce distortions at the level of 10%. 

Finally, we have compared our NLO QCD-corrccted results with existing calculations, 
which only take into account t/u-channel squarcd-diagram contributions. Working in this 
approximation, which renders the QCD corrections particularly simple, we found technical 
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agreement between our results and the existing calculations within statistical integration 
errors. We also found that, when typical VBF cuts are apphed, our full NLO QCD results 
agree with the ones in the t/u-channel approximation within fractions of a per cent. 

With the complete knowledge of NLO QCD and electroweak corrections, the theoret- 
ical uncertainty from missing higher-order effects should be of the order of 1-2% in total 
cross-section predictions for Higgs-boson masses in the range 100-200 GeV. For distribu- 
tions, the uncertainty will be larger in suppressed phase-space regions. The phenomeno- 
logical error of the parton distributions contributes a further 3.5% to the uncertainty, as 
reported in Ref. [17]. We thus conclude that the presented state-of-the-art results match 
the required precision for predictions at the LHC. 
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